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Abstract. Wc study the Wigner function for a quantum system witli a discrete, 
infinite dimensional Hilbert space, such as a spinless particle moving on a one 
dimensional infinite lattice. We discuss the peculiarities of this scenario and of the 
associated phase space construction, propose a meaningful definition of the Wigner 
function in this case, and characterize the set of pure states for which it is non-negative. 
We propose a measure of non-classicality for states in this system which is consistent 
with the continuum limit. The prescriptions introduced here are illustrated by applying 
them to localized and Gaussian states, and to their superpositions. 
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1. Introduction 

The formalism of Wigner functions and the formulation of quantum mechanics in phase 
space have been used since the early days of quantum physics. Originally motivated by 
the attempt to describe quantum effects in thermal ensembles, various quasi-probability 
distribution functions have been developed and applied to many different fields in 
quantum physics, as alternative formalisms which provide useful computational tools 
and facilitate physical insight into the quantum nature of states [U [2]. 

The first quasi-probability distribution function was introduced by Wigner in 
1932 [3j, in the context of statistical mechanics, to study quantum corrections to the 
thermodynamic equilibrium properties. In analogy to the classical situation, in which a 
state can be completely described in terms of its phase-space density, a quantum state 
can also be entirely characterized by its Wigner function, and the expectation values 
of all observables can be computed as a sum over the whole phase-space weighted by 
this function. In contrast to the classical case, in the quantum scenario probability 
distributions cannot be defined simultaneously over position and momentum. Thus the 
Wigner function is not a true probability distribution, as becomes apparent in the fact 
that it can adopt negative values. Instead, it can be interpreted as a quasi-probability 
distribution whose marginals reproduce the true probability distributions over single 
observables. Operators and dynamics can also be accommodated in the phase-space 
picture so that quantum mechanics can be entirely formulated in this framework. 

The Wigner function and other phase-space representations have found application 
to many different physical problems [1] , including quantum optics, statistical mechanics, 
hydrodynamics, nuclear theory and quantum field theory. In particular, in the field of 
quantum optics, the phase-space descriptions of quantum states have found extensive 
application. Specially interesting from the experimental point of view is the ability 
to reconstruct the Wigner function (and thus the quantum state) from measurements 
of the electromagnetic field quadratures, thus making it a very powerful tool for state 
tomography El C]- Remarkably, the fact that the Wigner function is not positive 
definite has itself a practical use, since the volume of its negative part can be used as a 
measure of non-classicality of the state [H]. 

In the last two decades, the interest in quantum information processing systems 
boosted the generalization of Wigner functions to quantum systems in finite dimensional 
Hilbert spaces. Although an early approach was pioneered by Stratonovich [H], who 
introduced a spherical, continuous phase space for a spin particle, more recent definitions 
target discrete phase spaces. The first such generalization was proposed by Wootters in 
1987 |10j for prime dimensional systems, and later generalized to any power of primes 
in |TT]. A different construction was followed in |T2l [131 E] which could cope with any 
dimension of the Hilbert space at the expense of enlarging the size of the phase space 
grid. The concept of negativity has also been studied in this context [I5l [161 E] 
connected to beyond-classical features of quantum algorithms. 

In this work we focus on a scenario not covered by the previous literature, namely 
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that of a quantum system with an infinite dimensional but discrete Hilbert space. The 
most immediate example would be a spinless particle moving on a one-dimensional 
lattice. We present a definition of the phase space and the Wigner function for this 
situation, which connects to both the features of the discrete constructions mentioned 
above and the proper continuum limit. 

The rest of the paper is organized as follows. In Sect. [2] we review the construction 
of the continuous and discrete Wigner functions used in the literature, and comment on 
the different features of each case. Sect. |3] presents our definition of the Wigner function 
for the infinite discrete quantum system, sketches the proof of its main properties, and 
shows that it reproduces the correct continuous limit. In Sect.|4]we introduce a measure 
of non-classicality of the states which can be computed from the Wigner function. Our 
definitions are illustrated in Sect. [5] with an explicit calculation of the Wigner function 
and the non-classicality for several examples. Finally, in Sect. [6] we discuss the utility 
of this definition, and how it can be applied to more general settings, for instance, to 
the case of a particle with spin moving on the lattice, or to several particles. 

2. Continuous and discrete Wigner functions 

The phase spaces of continuous and discrete quantum systems turn out to have striking 
differences. Defining a Wigner function for a discrete case, thus, requires more than a 
simple discretization of the continuum equations. There have been several prescriptions 
proposed for this kind of systems. In this section, after reviewing the main characteristics 
of the continuous Wigner function that a discrete version should respect, we summarize 
the different approaches that have been proposed, and their connection to the continuous 
definition. 

2.1. Continuous case 

For a quantum one-dimensional system, the Wigner function can be written [3] 

W,{x,p) = - / dy{x + y\p\x-y) e'^^P.^ (1) 

J-oo 

where p is the density matrix of the system, and |a;) represents the eigenbasis of the 
position operator, X [|] 

It is also possible to define the Wigner function axiomatically [3] . The fundamental 
properties that it must satisfy are usually formulated as follows. 

(i) Reality: the Wigner function is real. 

(ii) Projection: the integral of the Wigner function along any direction, in 
phase space, yields the probability distribution for the outcomes of measuring the 
observable aX + f3P, being P the momentum operator. In particular, then, the 

I We use natural units, such that h~l. 
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marginal distributions for position and momentum can be obtained, respectively, 
by / dpW{x,p) = P{x) and J dxW{x,p) = P(p).[|] 

(iii) Inner product: the inner product of two states, given by their density operators, pi 
and p2, can be computed from their Wigner functions as 

tr(pip2) = 27r J dxdpWi{x,p)W2{x,p). (2) 

The expectation value of any operator, O, can be also computed from its Wigner 
represent at ion , Wq {x,p). 

The Wigner function can also be constructed from the phase-point operators, defined 
for all points in the phase space as 

A{x,p) = -V{x,p)UV{x,p)\ (3) 

TT 

where V{x,p) are displacement operators and 11 is the parity reflection. The phase- 
point operators form a complete set, spanning all Hermitian operators. In particular, 
the Wigner function corresponds to coefficients of the density matrix in this basis, 

W {x, p) = ti {pA{x,p)) , (4) 

so that the full state can be reconstructed by 

p= dxdpW{x,p)A{x,p). (5) 



Properties equivalent to (|ij)-(|iiij) can be formulated for phase-point operators, 
leading to the same definition of the phase space. According to these properties the 
operators A{x,p) should be Hermitian and satisfy an orthogonality condition, and 
integrating A{x,p) along a line in phase space must yield a projector. 



2.2. Discrete finite case 

A valid generalization of the Wigner function to the case of a discrete Hilbert space 
involves generalizing the concept of phase space and the definition of phase space 
operators. Several approaches have been proposed in the literature for the case of a finite 
dimensional, periodic Hilbert space. Here we briefiy describe the two main alternatives, 
emphasizing their relation to the continuous case, and we establish the basis for our 
definitions (see [El E] for more comprehensive reviews). 

Wootters [TU] generalized the definition of the Wigner function to discrete periodic 
Hilbert spaces of prime dimension, N. In Wootter's original construction the phase space 
was a two dimensional N x N array, indexed by integers. Complete sets of parallel lines 
in this phase space, or striations, which are defined using arithmetics modulo A^, were 
associated to projective measurements. For more general cases, such as power of primes 
dimensions or composed systems with both discrete and continuous degrees of freedom, 

§ In the rest of the paper, unspecified integral (or sum) hmits will be understood as extending over 
the whole range of the integrated (summed) variable. 
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Figure 1. Graphical representation of the phase space for an infinite one-dimensional 
lattice. The momentum-like coordinate is continuous and periodic, k G [— 7r,7r[, and 
the position-like coordinate is discrete, labelled by integer values, m. 

the phase space could be constructed as a Cartesian product of the fundamental pieces. 
A related, more general construction, valid for systems whose dimension is an integer 
power of a prime number, was put forward in [TT]. In the general construction the 
discrete phase space has also size N x N, and was labelled by a finite field. To give 
a physical interpretation to the discrete phase space, each line is assigned to a pure 
quantum state. The set of all lines parallel to a given one corresponds to an orthogonal 
basis, and the distinct sets are mutually unbiased basis JT9\. This assignment of states to 
lines determines the particular definition of the Wigner function for the system, which is 
therefore not unique. Although closely connected to quantum information concepts and 
useful to describe systems of n qubits ISHIIIS], the lack of a unique physical interpretation 
and the restriction to dimensions which are powers of primes makes this approach less 
appropriate for the kind of system we want to describe. 

Leonhardt [13] introduced another definition, more closely connected to the 
continuous construction, in which the labels of the phase space axis could be connected 
to discrete position and momentum basis of the physical system. For the case of an odd 
dimensional system, the discretized version of the continuous definitions is enough to 
obtain a valid definition of the Wigner function. However, in the case of even dimensional 
Hilbert spaces, the naive discretization does not suffice to guarantee a Wigner function 
with the desired properties. Instead, half-odd labels had to be introduced between the 
integer points of the phase space axis, so that the size of the grid has to be increased to 
2N X 2N (see also [21] for a discussion). A similar approach was pursued in [14j, where 
the construction followed from the definition of discrete phase-point operators and was 
then applied to the analysis of quantum algorithms [22] . In [23] , this approach was also 
combined with Wootters' prescription to compose degrees of freedom and employed for 
the study of quantum teleportation. Our approach builds up on this construction, easily 
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connected to the physical interpretation of the continuous phase space. 
3. Definition in the infinite discrete lattice 

We consider here a single particle moving on a discrete one-dimensional lattice, with 
inter-site spacing a. We can define a discrete position basis, given by orthonormal 
states \n), with n E J.. Its Fourier transform defines then a quasi-momentum basis, 
I*?) ~ V^Xlra 'which can be restricted to the first Brillouin zone, q G [— f , f [• 

Unlike the discrete cases considered above, the Hilbert space of this system is not 
periodic and has infinite dimension. The continuum limit is recovered as a — )■ 0, and 
-^\n) — )■ |x = na). We may require that the Wigner function for this system, besides 
fulfilling the defining properties, also reproduces in that limit the usual one for a particle 
in one dimension, ([T]). 

In ^T3\ E] the problems in the direct discretization of Eq. ^ were connected to 
the fact that, in the case of even dimensions, such a definition does not generate enough 
independent operators. We may thus wonder whether the fact of having an infinite 
dimensional system is enough to solve this problem. However, the direct discretization 
of ([3]), in spite of producing an infinite number of operators, does not suffice either in 
this case to obtain a Wigner function that fulfills the desired properties in the case of 
interest. 

Here we follow closely the construction of [1A\, starting with a definition of the 
discrete phase space and the associated phase-point operators that then produce the 
Wigner function. 

A direct discretization of the phase-point operators leads to 

1 - 2k 

Adirect(m, k) = — f/2"HF( )e^2fcm^ (g) 

ZTT a 

where U is the discrete translation operator, shifting the lattice by one site, U"^\n) = 
\n + m), and V{q) is the continuous momentum translation, defined by its action on 
the momentum basis as V{q')\q) = \q + q')- Notice the hybrid character of the phase 
space in this case, with discrete and unbounded values of m and continuous, periodic 
k G [— TT, 7r[. It is easy to see that ^direct has periodicity vr in the momentum coordinate, 
^direct A; ± tt) = y4direct(^, k) . The Wigner function following from these phase-point 
operators does not fulfill the defining properties. In particular, summing over positions 
does not produce the correct marginal. Instead, 

5^iya.ect(n, A:) oc (^|p|^> + (^|p|^> . (7) 

n 

As a consequence, the resolution in k is not enough to retrieve all the information on 
the state from the Wigner function (see [131 EI] for a discussion of this effect). 

The problem does not appear if we integrate over the momentum coordinate, and 
is thus an effect arising purely from the discrete character of the position basis. It is 
not surprising, then, that the strategy in [H], consisting in doubling the number of 
points in the phase space, serves us to define also here an appropriate set of phase-point 
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operators. In our case, the doubling should only affect the position coordinate, and is 
equivalent to adopting the definition 

A(m,fc) = ^f/"W(-f )e*"^^ (8) 

In the position basis, the phase-point operators can then be written as A{m, k) = 
2^ Xln 1"^ ~ n)(n|e~**-^"~'^-"^, and the Wigner function for our system reads 

W{m, k) = tr [pA{m, k)] = — \"{n\p\m - r^)e-*(2'^-™)^ (9) 

n 

This corresponds to a phase space with the structure depicted in figure [T| where the 
m coordinate takes integer values, whereas k is continuous and periodic, taking values 
in [— 7r,7r[. Notice that (m, k) cannot be directly interpreted as discrete positions, and 
quasimomenta. Instead, they have to be understood as labels of the phase space points. 
To make the distinction clear, we reserve symbols m and k for the phase space and ra, 
q for the position and quasimomentum states of the lattice. 

Although the periodicity of A{m, k) is now 2-7?, it is important to notice that not 
all the operators are independent. Indeed, A(m, A; ± tt) = (— l)™A(m, fc), from which it 
follows 

W{m, k±n) = {-l)"'W{m, k). (10) 

It is easy to check that the definition ^ fulfills the main properties we require from 
a valid Wigner function. In particular it is real, as follows from the Hermiticity of ^ 



The inner product property (iii) is also easy to check, given operators A and B, 

27r ^ y dkWA{m,k)WB{m,k) = ii\^AI3y (11) 



m=— oo 



In a very similar way, we obtain the explicit expression of the density operator in 
terms of the Wigner function, 

/+7r 
dkW{m,k)A{m,k). (12) 



Due to the relation (jlOj), the orthogonality relation between phase-point operators 
adopts the following form, 

tr [A(mi, A;i)A(m2, k^)] = ^6^,m, [^(^1 - k^) + (-l)™^e(A;2)(5(A;i - k, + tc) 

+i-lr'Qi-k2)Sik^-k2-7^)], (13) 
where Q{k) is the Heavide step function. To obtain this relation we made use of 



^^ginfc _ 2nY2^S{k + 27rr), where the sum runs over all r G Z. Eq. (13) refiects 
the fact that operators associated to phase space points whose k coordinate is shifted 
by TT are not independent, but differ only in a phase. 

We may also compute the marginal distributions of (|9]), and obtain 

+00 

E W{m,k) = l{^Jp\^J, (14) 
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and 

r dkW{m,k) = Y,<2n{n\p\n). (15) 

The last equations make evident the distinction between the coordinates of the 
momentum space points, m ^ Z, k E [— 7r,7r[, and the position and quasimomentum 
bases, n, q. The k coordinate is adimensional and does not directly represent a 
momentum value, but is connected to q = k/a. The spatial label m in phase-space 
is only connected to a discrete position, s, for even values , m = 2s, while the odd 
values of m are analogous to the odd half- integer phase space grid points in [131 E]- 

Keeping these considerations in mind, we can take the continuum limit that 
transforms our discrete lattice into real space. This limit is attained by letting a — > 0, 
with na — )■ X G R. With this prescription, we can easily see that the continuum limit of 
Eq. ([9]) yields (up to a proportionality factor) the proper continuum Wigner function, 

W{m, k) -Wc{y = — , g = -) 

a-s>0 / 2 a 

1 /"^°° , ,ma , ,ma 



2tt y_oo ^ ^ 

as can be checked from the definition ([T| after a simple change of variable. Together 
with the discussion above, this result shows how the proper continuum limit is attained 
in the phase space coordinates. Indeed, as the spacing is decreased, m| — )■ x and 
^ — )■ g. The Wigner function is a quasi-probability distribution, and the physically 
meaningful quantities are given by integrals over the phase space. The measure of the 
integration must be modified according to this change of variables, so that we obtain 
the correspondence, for the integral over any region of the phase space. 



a-s>0 



dkW{m,k) — > dy dqW{y,q). (17) 



4. Non-classicality of states: negativity of the Wigner function 

The fact that the Wigner function is not positive definite over the phase space is 
interpreted as a quantum feature, since it follows from the incompatibility of quantum 
observables. This property has been applied to separate quantum states from classical 
ones. In the continuous case, it is known that the only pure states with non-negative 
Wigner function are Gaussian states [21]. The classification is not so clear for mixed 
states, where nevertheless some bounds are known for states with positive Wigner 
function [25|. From a quantitative point of view, the volume of the negative part of the 
Wigner function can be used as a measure of non-classicality [8] . More recently, it has 
been shown that the smallest distance to a state with positive Wigner function can also 
be used to measure the non-classicality of a state, without needing full tomography [26] . 

In the context of discrete systems the negativity of the Wigner function has also 
been explored, but the different prescriptions discussed in the previous section lead 
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to different conclusions. For the direct discretization, a discrete version of Hudson's 
theorem was proven by Gross [271 IT7j for the case of a Hilbert space with odd dimension. 
In that case the only pure states with non-negative Wigner functions are stabilizer states. 
For the class of discrete Wigner functions defined as in [11] a characterization was given 
in where the set of non-negative states was identified with the convex hull of 
stabilizer states when the Hilbert space dimension was small. Another scenario studied 
in the literature, which closely relates to our construction, is that of a pair of quantum 
variables, angular momentum and angle, and their associated phase space [281129] which 
has the same mixed discrete-continuous structure of figure [1} For that case, it has been 
shown that the only states with non-negative Wigner functions are those of well defined 
angular momentum. 

For the situation studied in this paper, a similar reasoning to that in [2H] leads to 
the conclusion that a pure state has a non negative Wigner function if and only if it is 
a state of well defined position in the lattice, i.e. the components of the state vector in 
the position basis are given by a delta function (see the Appendix for details). 

With our definition, however, phases (10) imply that any state with a non- vanishing 
Wigner function on some phase space point with odd-valued position-like coordinate, 
W{m = 2s + 1, k) ^ 0, will necessarily have a contribution of opposite sign at points 
(m, k ± tt). These signs are fundamental in order to ensure that the Wigner function 
reproduces the momentum and position probability distributions, but are not related 
to the quantumness of the different states. Therefore a naive calculation of the volume 
of the negative part of the function, i.e. applying the discrete version of the definition 
in [8] , will not be a valid measure of non-classicality, as it would result in a non- vanishing 
value even for states expected to be classical, such as the discrete version of Gaussian 
states. 

A similar phenomenon has been observed in different contexts. In the field of 
signal analysis, where Wigner functions have also been widely employed, the discrete 
time Wigner distribution shows similar features, which are related to aliasing [30], and 
various alternative definitions have been proposed to construct alias-free distributions, 
and to allow a reconstruction of the continuum time signal from a discrete sample. In 
the context of finite dimensional quantum systems, a proposal for a ghost free Wigner 
function was put forward in [31]. In all such cases, the negative values of the Wigner 
function respond to the very structure of the discretized phase space and not to the 
features of the state or the signal. We would thus like to define a new quantity which 
serves to estimate non-classicality of states in our system, and allows for a connection 
to the well-defined continuum limit. In particular, we expect that this non-classicality 
measure vanishes for all Gaussian states, so as to reproduce the well-known continuum 
limit, and that it does not include the spurious negative parts from the extended phase 
space. 

The definition (|9| leads to a discrete Wigner function which contains two images, 
one in each half of the momentum domain. According to (10), they have the same 
magnitude, but on odd position-like coordinates m, their sign is reverted, as can be 
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seen in figure [3] for tlie case of a pure Gaussian state. Although in the continuum hmit 
(|9| reduces to the original expression for the Wigner function, and the second, ghost 
image disappears, we would like to have a quantity that characterizes the non-classicality 
of discretized states. In particular, we require that the criterion is consistent with that 
for the analogous continuous states, in the cases when such exist. 

The so-called ghost image exhibits alternating signs between even and odd space- 



like coordinates (see, for instance figure 4(a)), while the regular image is smooth. We 



would thus like to use as a measure of the non-classicality the negativity restricted to the 
regular image. However its position in the phase space plane is not fixed, but changes 
with momentum shifts. A momentum displacement, go; translates into a displacement 
qoa in the k coordinate. As the lattice spacing vanishes the regular image lies on the 
central region of the phase space, while the ghost image is pushed towards the edge, 
which in the continuum is mapped to infinity. 

Instead of trying to locate the regular image, so as to restrict the sum to the 
corresponding phase-space region, we may apply a filter that eliminates the spurious 
sign oscillations from odd values of m, and effectively produces two copies of the regular 
image. We thus define the following quantity, 

+00 

r]{p)= J2 / dk[\W^'\m,k)\-W^'\m,k)], (18) 

J — 7r 

m=— 00 

where ly^*) is the result of filtering out the sign oscillations for odd m. If the filtering 
is perfect, in the continuum limit ri{p) will yield twice the negativity of the Wigner 
function as defined in [8]. 

Different filtering operations can be tried to this aim. In particular, we propose to 
use a sign-averaged Wigner function, defined by 

lyWf^ M = / ^^""'^^ (19) 
^ ~ \ xim,k)\Wim,k)\ modd, ^ ' 

with x('^) k) = sign[2 sign(l^(m— 1, k)) + sign(W{m, k))+2 sign(W{m + l, k))], i.e. the 
even components are unchanged, and the sign of the odd ones is corrected according to 
a majority criterion that takes into account the sign of the two closest neighboring even 
points [II] This produces approximately two copies of the regular image (see figure 4(b) ), 



so that ri{p) is equivalent to twice the negativity restricted to the half space where this 
image is supported. 

5. Particular cases 

To illustrate the definitions introduced in the previous sections, we explicitly compute 
here the Wigner functions and negativities for several pure states. 

II The factor 2 takes care of the situation when one of the adjacent even points has vanishing W{m, k). 
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Figure 2. Wigner function for a state localized at the origin uq = 0, assuming a = 1. 
Notice that a plane projection of the phase space is represented, although it is periodic 
in k, and the edges k = ±7r have to be identified. 



5.1. Localized state 

We may consider the simplest case, in which the position of the particle is well defined, 
l^is) = I'^o), so that, in position basis, (nl'^s) = ^nno- The Wigner function can then be 
computed exactly, 

Ws{m,k) = ^6m2no- (20) 

Zn 

This function, represented in figure [2| is non-negative everywhere, so that r](5) = 0. 

5.2. Gaussian state 

In the case of continuous degrees of freedom, Gaussian states play a fundamental 
role. In particular, pure Gaussian states are the only pure states with non negative 
Wigner function [23]. It makes then sense to consider the discretization of a state 

^ _ (x — xq)^ , _ (n — nQ)2 

\E'(x) = — ^-T-e 2<t2 g^m^^ namely \^g) = ;^S„e 2*^ e*''°""|n), for uq E Z being 

a = a/a the width measured in units of the lattice spacing. The correct normalization 

in the discrete case, A^^ = ^ = 6*3(0, e~5^), is expressed in terms of the 

Jacobi theta function, defined as 9^{z.,q) = ^^g^^e^*^" for complex arguments q, z, 
with |g| < 1 [32j. 

The Wigner function for this state can also be computed exactly. 



1 



WG{m,k) = _e*(fc-9oa)mg ^ ^" Y (21) 

27r ^3(0, e-^) 

and shown in figure [3] for the particular case a = 2, uq = 0, qq = 0. The figure shows 
clearly the regular image, centered around k = 0, and the ghost image, exhibiting the 
sign oscillations on odd sites. If we consider instead a displaced Gaussian, with go 7^ 0, 



the whole figure is correspondingly shifted in momentum space, as shown in figure 4(a 
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Figure 3. Wigner function for a discretized Gaussian state with a — 2, qq = and 
uq = 0, taking a = 1. 




Figure 4. Wigner function (left) for a discretized Gaussian state identical to that in 
fLgure|3]with a displacement in momentum qqg = 7r/3. On the right, the sign-averaged 
Wigner function (19). 



To illustrate the meaning of the sign-averaged function defined in (19), we also plot it 



in figure 4(b) for this state. Obviously, ri{'^G) = for any pure Gaussian state. 



5. 3. Superposition of deltas 

The Gaussian case has vanishing negativity, as expected from its correspondence in the 
continuum limit. It actually includes the case of a localized state, too, which can be 
interpreted as a Gaussian in the limit of a vanishing width, a. Superpositions of such 
states will instead have more quantum features. 

We may in particular consider an arbitrary superposition of two localized states, 
such as \'^25) = ^il\^\2 ('^""1 + '^^nnz) I""-), ^OT any ni 7^ n2 G Z and a G C. The 
corresponding Wigner function can be easily calculated, 

zvrl 1 + ar j 
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Figure 5. Wigner function for the superposition of two deltas located at ni = — ^2 = 
4, assuming lattice spacing a — 1. The function vanishes everywhere except on three 
isolated strips, colored purple in the figure. 



+2\a\6m,ni+n2 cos[An k + 4>]} , (22) 

where is the phase of the complex coefficient a, and An = n2 — ni. In this case, the 
Wigner function vanishes everywhere except for three particular values of the space-like 
phase space coordinate, namely m = 2ni, 2^2, ni + n2. Figure [5] shows the particular 
case of ni = —n2 = 4, a = 1. 

It is easy to see that W2s\m, k) = W2s{'m, k), since none of the terms changes sign 
under (19). Indeed, the ffist two terms are non- vanishing only on even values of m, 
while the last term can be supported on odd m if ni and n2 have different parity, but 
in that case, W2s{fn ± 1, /c) = 0. Therefore, we can analytically compute the quantity 
^ as 

/TT 
[\W{m, k)\ - W{m, k)]dk 





a 




7r(l + 


a 





/TT 
dk [\ cos(An/c + 0)| — cos(An/c + 



4 


a 




7r(l + 1 


a 





(23) 



independent of the separation between the localized states. An, and reaching its 
maximum value, 77max(*2i) = S/vr, for |q;| = 1. 



5.4- Superposition of Gaussian states 

Another family of states for which the Wigner function defined above can be computed 
analytically is that of superpositions of pure Gaussian states. We may consider an 
arbitrary superposition of two discretized pure Gaussian states. 




Wigner function for a particle in an infinite lattice 



14 



for arbitrary ni^2 £ Z, qi^2 G [— vr/a, 7r/a[ and a G C. For such state, the Wigner function 
can be expressed as a sum 

W2G = Wi + \a\^W2 + aWu + a*W2u (25) 

where Wi and W2 are (up to the normahzation factor) equivalent to the Wigner function 
of a single Gaussian (21), while 14^12 and 14^21 contain the crossed terms, 

1 



Wuim, k) 



(m-n2) 



X 6'q A; — a 



qi+q2 



+ i 



2^2 2^2 



'^l+°^2 



(26) 



and W21 = Wi2{l ^ 2). 

In the symmetric case, a = 1, rii = — n2 = no, ai 
expression adopts the compact form 



a2 = cr, gi = g2 = 0, the above 



W2G{m, k) 



"0 mno 
e s^cosh— — h cos(2A;noj 



^2 



m 



(27) 



with A/"^ = 2(1 + e "o/'^^)03(O, e ^^'^'^). Using the properties of the 63 function, we can 
further simplify the expression, so that, for even m = 2s, 



W2Gi2s,k) 

and for odd m = 2s + 1, 



.2 



TT. 



-Af- 



j_ . 
^2 



e g^cosh ^ ° + cos(2A;rio) 



(2J 



iy2G(2s + l,A;) 



(s+l/2)^ 



7rA/'2 



2^2 



4 (2s + lH^ , 
<^2 cosh — h cos(//enoj 

(j2 



(29) 



In the limit a — )■ 0, (28) results in the expression for the superposition of two 
localized states discussed in the previous section, while (29) vanishes. 

Figure 6(a) shows the full Wigner function for the particular case no = 6, a = 1.5. 
The central part, around k = 0, corresponds to the regular image, showing the usual 
Gaussian peaks and a central interference region. This survives in the continuum limit, 
giving rise to a genuine negativity. The ghost image in this case lives on the half phase- 
space with larger momenta, and exhibits the characteristic sign oscillation when moving 
along the space-like axis. The sign average defined in (19) transforms this image in 
a copy of the genuine one, as shown in figure 6(b) , so that ?7(\I'2g) will be twice the 
negativity of the regular image. 

Although there is no closed analytical expression for the non-classicality ri{'$2G), 
even in the simplest case discussed above, we can compute it numerically, as shown in 
figure 7(a) for the symmetric superposition of two Gaussian states of the same width, 
centered at ±nQ and with momentum displacements, gi = and q2 = qo- As shown 
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Figure 6. Wigner function (left) for the symmetric superposition of two Gaussian 
states with center in ±no, for uq — 6, and width a = 1.5, assuming lattice spacing 
a = 1. The right panel shows the sign-averag ed for comparison. 




Figure 7. Negativity of a superposition of two discretized Gaussian states of the same 
width. The left plot shows the case a = 1.2, as a function of their half- distance, no, 
and the relative momentum displacement, qq, taking a = I. On the right, we show the 
asymptotic value reached for varying width, a. 



in the plot, rj vanishes only for no = and go = 0, vr, when the situation reduces to a 
single Gaussian. For small distances, 2no, the value of rj depends on no and Qq, while for 
larger separations it becomes less sensitive to qo, and soon enough it reaches its maximal 
value, and stays constant. As shown in figure 7(b), this asymptotic value is sensitive to 
the Gaussian width only when the latter is comparable to the lattice spacing. When a 
is large enough, instead, the asymptotic negativity is constant. In the limit a — 0, on 
the other hand, the negativity for a superposition of two deltas is recovered. 



6. Discussion 



We have extended the formalism of the Wigner function to the case of a quantum 
system with a discrete, infinite dimensional Hilbert space. For instance, this would be 
the case for a spinless particle moving on a one dimensional lattice. The prescription 
presented here appears to be the natural one for this problem, as it satisfies the defining 
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mathematical properties of the phase-space representation and recovers the correct 
continuum hmit for vanishing lattice spacing. 

The quantification of non-classicality, as signaled by the negative part of the Wigner 
function in the case of continuous degrees of freedom, has to be redefined in this case 
to exclude negative contributions due to the structure of the discrete phase space itself. 
We have proposed a negativity measure for this case, and have illustrated it with the 
explicit results for localized and Gaussian states, and for superpositions of each. Our 
results support the meaningfulness of this measure to characterize the states of a particle 
on a one-dimensional lattice. 

As for other cases, in which the phase-space formalism can also account for the 
dynamics of the system, it would be possible to formulate the evolution of such system 
fully in terms of its Wigner function, and to use the proposed measure, //(p), to classify 
quantumness in evolving states. Although the examples presented in this paper are 
focused on pure states, the same concepts apply also to mixed states. 

An interesting extension of this work is combining the phase space introduced here 
with additional degrees of freedom, such as internal ones for the particle, or to extend 
it to the case of several particles or dimensions. Wootter's prescription [TU] to construct 
composite phase spaces by combining the phase-space point operators of different degrees 
of freedom via their tensor product can be applied in this case. 
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Appendix A. Pure states with positive Wigner function 

Analogous to the result in [28j for a conjugate pair of angle and angular momentum 
variables, with the present definition the Wigner function of a pure state is non negative 
if and only if it is an eigenstate of the discrete position operator, i.e. {n\^) = 6nno- The 
first part of the theorem is trivial, since the Wigner function of a localized state (20) is 
non negative. 

To show the converse, let us assume a pure state with non-negative Wigner function. 



W{m,k) > 0, Vm G J., k G [— vr, 7r[. From (10) it follows that the Wigner function can 



only be non-vanishing on points of the phase space with even space-like coordinate, 
m = 2n, 

W{2n + l,k) = WneZ. (A.l) 

The rest of the demonstration follows closely that in [28], and we sketch it here 
only for completeness, with the proper modifications to match the definition in (|9|). 



Wigner function for a particle in an infinite lattice 



17 



The proof relies on the following two lemmas, proven in [28|, for complex periodic 
functions and their (discrete) Fourier transform. 

(i) Let g{q) be a continuous, complex, 27r— periodic function. If its Fourier transform 
is non-negative, then the integration kernel g{q — q') is non-negative. 

(ii) Given a function / : Z — )■ C, if its inverse Fourier transform has constant modulus, 
then Y.n&z fin)f*{n + m) = 0, Vm ^ 0. 

It is easy to see that, for a pure state, the Wigner function can be written as 

W{m, k) = ^ r''^ dge^«"^>(f + q)r{l - q) (A.2) 

J-n/a 

where ip{k) = (klip) are the components of the state in the quasi-momentum basis. It is 
thus the Fourier transform of the function g{q) = ^J^ -?/^((fc + q)/a)ip*{{k — q)/a). From 
lemma i| f^^ d<l'x*{Q)9{Q~ Q')x{q') ^ for any x- Iii particular, requiring the inequality 
for all functions x{<l) = aiS27z{q - Ci) + a262n{q - C2), where 62n{q) = Z]rez^(9' " 2r7r), 
0-1,2 £ C, ci,2 € IR, implies 

IV^(g)r> IV^(g + A)||V^(g-A)|, vAgR, (A.3) 

This requires that li'iq)] is constant, so that also \g{q)\ must be constant. Applying 
now lemma ^ 

^Vr(m,A;)iy(m + j,A;) = ^ 0. (A.4) 
j 

So that, for a given value of k, there can at most be a single space-like component, 
mo{k), for which the Wigner function does not vanish. Combining this with (A.l), we 
obtain that such component must be even, mo (A;) = 2nQ{k), so that, using normalization, 

W{m, k) = ^Sm,2no{k)- 

It only remains to be shown that this component is the same for all values of k. 
This can be seen, as in [28], by making use of the expression for the Wigner function 
for a product in terms of individual Wigner functions. 



1 

{m,k) = — / dki / dk2Wo^{m + mi,k + ki) (A.5) 

271" J^n J-n 



mi,m2 



W,,{m + m2,k + k2)e''^'^^^'-'^'^^\ 

In particular, taking pi = p2 = p, the pure state we are considering, for which 
W{m,k) = ^Sm,2no(k), and looking at the (real) component for m = 2no(0) = 2nQ, 
k = 0, 

dki / dk2 cos {2k2[nQ{ki) — Uq] — 2ki[nQ{k2) — no]) . (A. 6) 

-TT d —TV 

To fulfill this equality, the argument of the cosine has to be an integer multiple of 27r 
for all values of ki^2, which is only possible if nQ{k) = uq V/c, and thus 

iy(m,A;) = -^5„,2no- (A.7) 

ZTT 

Using ( 12 ) it is easy to show that the pure state corresponding to this Wigner function 
is \^s) = l^o)- □ 
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